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1. Introduction 

A diffusion-like description of radiative transfer in 
clouds and the free atmosphere is often employed. 
The two stream model is probably the best known 
example of such a descript ion. 6 ’ 12,9 The main idea 
behind the approach is that only the first few mo- 
ments of radiance are needed to describe the radia- 
tive field correctly. Integration smooths details of 
the angular distribution of specific intensity and it 
is assumed that the closure parameters of the theory 
(diffusivity factors) are only weakly dependent on the 
distribution. In this paper we investigate the diffusiv- 
ity factors using the results obtained from both Stra- 
tocumulus and Cirrus phases of FIRE experiment. A 
new theoretical framework is described in which two 
(upwards and downwards) diffusivity factors cure used 
and a detailed multi-stream model is used to provide 
further insight about both the diffusivity factors and 
their dependence on scattering properties of clouds. 

2. Diffusivity factors 


of specific intensity I(r } fx). Vector S' is related to 
the source term. The sign convention is as follows: 
the instrument located on a airplane and facing the 
Earth’s surface will measure the upward flux of radia- 
tion, here indicated by the + sign, the corresponding 
6 values are between (tt/2, 7t) and fx £ (0,-1). The 
instrument facing skywards will measure downward 
flux, here indicated by the — sign, the corresponding 
9 values are between (0,tt/2) and fx £ (1, 0). 

Equation (2.1) is a direct analog to the Fickian 
diffusion law where the flux quantity is related to 
the gradient of a scalar quantity (scalar irradiance) 
through the diffusivity matrix D. In oceanography 11 
and atmospheric science 2 another diffusivity matrix 
is often employed 

U = DF (2.4) 

where D is the 2x2 matrix 



and the diffusivities D and D+ are defined as 


There are many diffusion- like approximations in ra- 
diative transfer theory. 10 The most intuitive being 


U-(t) 

F~(t) 


and 


D + = 


F+(r)’ 


( 2 . 6 ) 


- d\J 

F = D— +S' 
dr 

where flux and scalar irradiance vectors 

F 


( 2 . 1 ) 


-(£)• u =(£ ! ) < 2 - 2 > 

are defined by the following hemispheric averages 

F + = J nI(T,n)dn F~ = J nl(r,/j.)di4,2.3a) 
U + = J I(t, n)dn U~ = J I(T,n)dn (2.3b) 


It can be seen that (2.6) involves the first two mo- 
ments of angular radiance dependence. These quanti- 
ties are used to close the set of hierarchy of equations 
for various moments of radiance. The resulting set of 
equations is called the two-stream approximation. 

3. Measured diffusivities from the 
CSU bugeye 

The Marine Stratocumulus Intensive Field Obser- 
vations (MSIFO) of the First ISCCP(International 
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Figure 1: CSU bugeye instrument 



Wnvetength [^im] 

Figure 2: Spectral response of CSU Bugeyes 


Satellite Cloud Climate Program) Regional Experi- 
ment (FIRE) was conducted off the coast of Cali- 
fornia at and in the vicinity of San Nicolas Island 
in July 1987. On San Nicolas Island, a tethered 
balloon with radiometric instrumentation from Col- 
orado State University (CSU) probed the marine 
stratocumlus. 7 The shortwave radiometric instru- 
mentation measured both the downwelling and up- 
welling irradiances, and consisted of four Epply pyra- 
nometers (measuring wavelengths from 0.3 to 2.8 jim 
and from 0,7 to 2.8 /im), and two CSU bugeyes. The 
bugeye measurements will be used to derive the ra- 
diative diffusivities. 

The Bugeye, the CSU Multidirectional Photodiode 
Radiometer 4 measures the angular distribution of the 
radiance field. It consists of a hemispherical array of 
thirteen silicon photodiodes and associated electrical 
circuitry mounted on an aluminum housing (Fig. 1). 
The upward looking bugeye had diodes with a 10° 
field of view. Each diode of the downward looking 
bugeye had a 50° field of view. The spectral range of 
both bugeyes is from 0.36 to 1.10 /rm but the diodes 
of the downward looking bugeye were covered with a 
blue tinted Schott glass filters. Figure 2 shows the 
spectral response of the downward looking bugeye 
(solid line) and the upward looking bugeye (dashed 
line). The peak sensitivity of the downward looking 
bugeye is at 0.40 /rm while the upward looking bugeye 
has a peak sensitivity at about 0.93 jim. The bug- 
eye voltages are actually irradiances seen within each 
bugeye diode’s field of view. These voltages were nor- 
malised to the sensitivity of the first diode after being 
corrected for a zero offset. A field of view correction 
was then applied to each diode voltage to provide the 
proper steradian weighting on the hemisphere. A cor- 
rection for the pitch and roll of the platform was also 
accounted for. From the 13 diode measurements, the 


diffusivity factor was deduced as 

P _ I Vj>l<(hgmi)/i4(diode) 

£i*i^>Mhemi)M(diode)cos0 4 ( ' ’ 

where is the voltage of diode i, A(diode) is the 
steradian field of view of each diode, A<(hemi) is the 
steradian area that diode i voltage represents, and 
is the zenith/nadir angle (accounting for the pitch 
and roll of the platform) of diode t. Comparison 
of (3.1) to (2.6) and (2.3) emphasizes the approxi- 
mations used to integrate over angle. To test the 
approximation, the irradiances obtained from sum- 
mation of the diode measurements (the denomina- 
tor of (3.1)), were compared to the shortwave irra- 
diances obtained directly from the Epply pyranome- 
ter. An excellent relationship (not shown) between 
these measurements was obtained despite the differ- 
ent spectral characteristics of the bugeyes and pyra- 
nometers. The diffusivities for the downwelling ra- 
diation (D~) and for the upwelling radiation (D + ) 
are shown in Figure 3 for a Sc cloud smapled on 
the morning of Julian Day 189 (8 July 1987). The 
diffusivity values average about 1.62 with D + being 
slightly larger than D~ . Cloud base is at about 970 
mb and cloud top is at about 930 mb. As the bug- 
eye entered the base of the cloud the value of D + 
rapidly increases until D + assumes a in-cloud profile 
and then decreases through out the cloud. D~ in- 
creases with height through the cloud to cloud top 
where it has has the same value as £> + . 

We also performed some preliminary calculations of 
diffusivities for the downward looking (and only) air- 
craft bugeye from FIRE Cirrus, Oct. 28, 1986. The 
results for several passes through two cirrus samples 
give results of D + « 1.7. 
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Figure 3: DifFusivity coefficient. Marine stratocumu- 
lus . FIRE 8 July 1987 


4. Model results 


A. Diffusivity factors from the de- 
tailed radiative transfer code 

To demonstrate further the structure and depen- 
dency of diffusivity factors under various condi- 
tions, we have employed a comprehensive radiation 
model 18 which includes scattering and absorption by 
molecules, as well as by particles. The midlatitude 
summer® atmospheric profile is employed in which a 
cloud layer of 1 km thick is located with a base 1 
km above the surface. The cloud layer is assumed 
to be homogeneous with a mean liquid water con- 
tent of 0.2gm~ 8 and a mean effective radius of 10/im. 
Three spectral bands of 0.52 — 0.57/xm, 1.28 — 1.53/11X1, 
and 2.38 — 2.91/un are selected, respectively to repre- 
sent the cases of no, moderate and strong absorption 
by water droplets. Three solar zenith angles of 0°, 
30° and 75° are also chosen to investigate the angu- 
lar dependence of diffusivity factor. The results of 
these calculations are shown in Fig. 4. These cal- 
culations clearly demonstrate the dependence of 
and D~ on absorption strength, solar zenith angle, 
and optical depth. The diffusion domain in which 
D + (r,/io) = D+ and D~(r t fio) = D~ is also appar- 
ent for large r. The results also indicate that the 
diffusivity is higly variable for clouds of r < 1. 


B. Two-stream approximation 

The azimuthally and hemisphericaly averaged 
monochromatic radiative transfer equation for diffuse 
radiation, 7(r, /i), in a plane- parallel, horizontally ho- 
mogeneous medium which scatters, emits and absorbs 
can be written as 

= u+ i r ) - 11,0 J 0 

-Wo [ -n)dn- S + (t)(4.1) 

Jo 

dF J^ m = ~ U ~( T ) P(n)I{r,-n)dn 

+u>o [ <p{v)I(T,ii)dn + S~(t)( 4.2) 

Jo 

where the backseat ter fraction is 

m = \ J P{th~^)dn' (4.3) 

and the forward scattering fraction 

<p(n) = l-P(n) = ~ J P(n',n)dn' (4.4) 

Equation (4. 1-4.2) doesn’t involve any approxima- 
tions. Unfortunately it is a system of 2 equations and 
6 unknowns: F* + , F' - , U“ and 2 independent 

integrals involving scattering fraction. To overcome 
this difficulty all two-stream models introduce some 
kind of closure hypothesis. 9 * 1 The particular choice 
below follows the discussion given in Preisendorfer 11 
and Buglia. 1 To close the system of equation we need 
the relationships between F, U (eq. 2.1 or 2.1, for in- 
stance), and the integrals involving backscatter. For 
this purpose we introduce diffusivities as defined by 
(2.6) and assume that the back and forward scatter 
fractions are independent of angle. It follows that 

dF 

= AF + S (4.5) 

dr 

where A is the 2x2 matrix 

a=(;( :£) <«> 

and 


it = D + (1 - G) 0 tp) , 7 i=D (l-w 0 <fi) (4.7a) 

^t+ = D + woP, — D~C>oP (4.7b) 


and 




(4.8) 
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Figure 4: Diffusivity factors from detailed radiative transfer model 









Matrix A depends on local optical properties of the 
atmosphere (such as the single scattering albedo u> 0 
and the asymmetry parameter g). With the typi- 
cal assumption that D+ = D~ — D> i.e. that the 
properties of the homogeneous layer are radiatively 
isotropic, we write scalars 71 and 72 as 

7! = D(1 - aw) (4.9a) 

72 = Du>q(3 (4.9b) 

For the case = D~ , the matrix A exhibits spe- 
cific structure and is defined by two scalars. It is this 
fact which allows us to introduce only two physical 
parameters defining the homogeneous and isotropic 
medium: transmittance t and reflectance r. For 

the anisotropic (but homogeneous) layer, for which 
£>+ ^ D~ t the matrix structure (4.0) is defined by 
four scalar entries and two transmittances and two 
reflectances are needed to define the system. The ap- 
parent 11 anisotropy of the medium is forced by dif- 
ferences in the distributions of the intensity incident 
at cloud top and base. The solution of (4.5) may be 
written in the form 5 

F (r) = PF (r t ) (4.10) 

where P is a 2 x 2 fundamental solution matrix. The 
propagator for a homogeneous atmosphere, i.e. for 
constant A , is 

F (r) = e A ( r-T ‘)F (r t ) (4.11) 

for a layer of thickness r — r t . For the case of a two- 
stream model and we can represent the propagator 
as 

Pmp= a7 

where 1 is the 2x2 identity matrix. The eigenvalues 
Ai and A 2 are roots of the characteristic equation 

A 2 -pA + <? = 0 (4.13) 

and 

p = tr A = 7^ - 7f (4.14a) 

q = det A = y~ — 7 + (4. 14a) 

It can be shown that the reflectances and transmit- 
tances are directly related to the propagator 

/ r+ t+ \ _ _1_ / ~Pi2 1 

\t“ r~ ) Pn \PiiP 32 - P21P12 P21 

where 

Pn = y [e A * T (7! + - A,) - e*’ T (7i + - *0] (4.1«a) 



i— [e x ' T (A - Ajl) - e x ’ r (A - X 1 1)] 
— A 3 

(4.12) 
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Figure 6: Reflectivities from the four parameter 

two-stream model 

Pn = -y 72 - (e x ' r - e x ’ T ) (4.16b) 

P2 i = y7 2 + (e X,T -e A,r ) (4.16c) 

P2 , = - j [« X,r (7f + A 2 ) - c X ’ T (7f + A x )] (4.16d) 

and l = Ai — A 3 . Figure 5 presents results for two 
sets of diffusivities. First = 2.2, D~ ~ 2 which 
corresponds to the no absorption case (compare Fig. 
4). The second set D+ = 2.8, D~ = 1.2 is more 
typical of the absorption case. The results indicate 
that the medium becomes anisotropic (r+ ^ r”) and 
that reflectances and transmittances of the cloud are 
sensitive to changes in the diffusivity matrix. 

5, Summary 

The diffusivity factor has been studied in the con- 
text of two-stream approximation of the radiative 
transfer equation using data obtained from both Stra- 
tocumlus and Cirrus phases of FIRE project. De- 
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tailed radiative transfer calculations have been em- 
ployed. The theoretical framework of the two-stream 
model for homogeneous but anisotropic (anisotropy 
being forced by the boundary conditions) atmo- 
spheric layer is described. It is shown that the to- 
tal extinction matrix is defined by four components. 
Preliminary results indicate that the diffusion coef- 
ficients depend on the specific spectral region un- 
der consideration, absorption strength, optical thick- 
ness and solar zenith angles. In the strong absorp- 
tion case, upward and downward diffusivities separate 
and the two-stream model predicts large differences 
in upward and downward reflectances and transmit- 
tances for such a case. Model results indicate the 
existence of a diffuse region in which the diffusion 
matrix doesn’t vary with respect to solar zenith an- 
gle and/or increased optical thickness. Experimental 
data indicate the variability with the sun’s zenith an- 
gle of the downward diffusivity in the free air overly- 
ing the cloud. In the cloud layer the upward diffusiv- 
ity is larger thanjjie downward diffusivity which is in 
agreement with the strong absorption case obtained 
from the numerical model. Downward diffusivity is 
fairly constant throughout the cloud and below the 
cloud layer. The upward diffusivity adjusts rapidly 
to the existence of the cloud layer. Model results 
indicate that a ‘skin’ layer in which rapid changes 
in diffusivity occur is of the order of r = 1. Thus 
constant- D models may not be suitable for thin cir- 
rus clouds. 
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